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ABSTRACT 

UM 625, previously identified as a narrow-line active galactic nucleus (AGN), actually exhibits broad Ha and 
H/3 lines whose width and luminosity indicate a low black hole mass of 1.6 x 10 6 M . We present a detailed 
multi wavelength study of the nuclear and host galaxy properties of UM 625. Analysis of Chandra and XMM- 
Newton observations suggests that this system contains a heavily absorbed and intrinsically X-ray weak (a ox = 
-1.72) nucleus. Although not strong enough to qualify as radio-loud, UM 625 does belong to a minority of 
low-mass AGNs detected in the radio. The broad-band spectral energy distribution constrains the bolometric 
luminosity to Lb i ~ (0.5-3) x 10 43 erg s _1 and Lboi/^Edd ~ 0-02-0.15. A comprehensive analysis of Sloan 
Digital Sky Survey and Hubble Space Telescope images shows that UM 625 is a nearly face-on SO galaxy with a 
prominent, relatively blue pseudobulge (Sersic index n = 1.60) that accounts for ~60% of the total light in the R 
band. The extended disk is featureless, but the central ~ 150-400 pc contains a conspicuous semi-ring of bright, 
blue star- forming knots, whose integrated ultraviolet luminosity suggests a star formation rate of ~0.3 M yr _1 . 
The mass of the central black hole roughly agrees with the value predicted from its bulge velocity dispersion but 
is significantly lower than that expected from its bulge luminosity. 

Subject headings: galaxies: active — galaxies: individual (UM 625) — galaxies: nuclei — galaxies: Seyfert — 
X-rays: galaxies 



1. INTRODUCTION 

Supermassive black holes (BHs), with masses in the range of 
10 6 - 10 10 M as measured via stellar and gas kinematics, have 
been convincingly inferred to be present in the centers of nearby 
inactive massive galaxies and are generally believed to reside in 
all galaxies with a spheroidal stellar component (see Kormendy 
& Ho 2013 for a review). Moreover, there are tight relations be- 
tween the BH mass (Mbh) and the properties of the spheroidal 
component (namely, ellipticals and the bulges of disk galaxies), 
including stellar velocity dispersion (a*; Ferrarese & Merritt 
2000, Gebhardt et al. 2000a), luminosity (Lbuige; Kormendy & 
Richstone 1995; Magorrian et al. 1998; Marconi & Hunt 2003), 
and mass (Mb u i ge ; Haring & Rix 2004). 

However, the situation is far from clear in the low-mass 
regime (M B h ^ 10 6 M ) because such BHs are largely be- 
yond the reach of current capabilities for direct dynamical mea- 
surement. A practical approach is to search for them in low- 
luminosity type 1 active galactic nuclei (AGNs) with mass esti- 
mated from their broad-line width and luminosity using empir- 
ical virial relationships (e.g., Gebhardt et al. 2000b; Kaspi et 
al. 2000). This technique has yielded a sample of ~ 200-300 
candidate BHs with masses between 10 5 to 10 6 M (Greene & 
Ho 2004, 2007b; Dong et al. 2012b), including a couple be- 
low 10 5 M , in the regime of so-called intermediate-mass BHs 
(IMBH, 10 3 " 6 M ; Filippenko & Ho 2003; Barth et al. 2004; 
Dong et al. 2007; see reviews in Ho 2008; Kormendy & Ho 
2013). The host galaxies of low-mass BHs thus found appear 
very different from their supermassive counterparts. The best 
nearby example, NGC 4395, is a dwarf Sdm galaxy without a 
bulge at all (Filippenko & Ho 2003), while the second proto- 
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type, POX 52, is instead a spheroidal or dwarf elliptical galaxy 
(Barth et al. 2004; Thornton et al. 2008). According to the 
studies to date, the M B h-o"* relation of local inactive massive 
galaxies appears to roughly extend to the low-mass end (Barth 
et al. 2005; Greene & Ho 2006; Xiao et al. 201 1), but the M BH - 
^buige relation does not. Photometric decomposition by Greene 
et al. (2008) and Jiang et al. (201 la) indicates that the majority 
of the host galaxies with disks are likely to contain pseudob- 
ulges; the rest resemble spheroidals according to their position 
on the fundamental plane. Very few live in classical bulges. 
Moreover, the M B H-^buige relation flattens out at the low-mass 
end, and on average the bulge luminosity is larger by 1-2 index 
at fixed M BH (Greene et al. 2008; Jiang et al. 201 lb). 

The low BH masses and high accretion rates that generally 
typify the low-mass AGN sample also provide a unique oppor- 
tunity to probe accretion processes in an under-explored regime 
of parameter space. Especially interesting is the broad-band 
spectral energy distribution (SEDs) of these systems. Unlike 
quasars, which have luminosities that far exceed that of their 
host galaxies, the SEDs of Seyferts, particularly those with 
low-mass BHs, can be strongly contaminated by emission from 
their hosts. Observations with high angular resolution are ab- 
solutely essential to decouple the AGN from the host to con- 
struct proper nuclear SEDs. In spite of these challenges, stud- 
ies so far already suggest that low-mass AGNs may possess 
unusual multiband properties, including for the tendency to be 
very radio-quiet (Greene et al. 2006; Greene & Ho 2007b) 
and relatively X-ray bright (Greene & Ho 2007a; Miniutti et 
al. 2009; Desroches et al. 2009; Dong et al. 2012a). 

The most extensive multiwavelength studies on low-mass 
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Fig. 1. — The SDSS spectrum of UM 625 (black), together with its continuum decomposition and line fitting. (Left) The model for the starlight, featureless 
continuum, and the sum of the two are plotted as green, blue, and red lines, respectively. The residual spectrum (observed - model) is shown on the bottom, offset 
downward for clarity. The inset highlights the H/3 and [O III] region and the clear detection of broad H/3 in the residual spectrum. (Right) Decomposition of the 
Hck+[N II] region. The upper panel shows the original data (thick black line), the fitted narrow lines (green), the fitted broad Ha component (blue), and the sum of 
all the fitted components (red). The lower panel shows the residuals of the fit. 



BHs thus far have concentrated, unsurprisingly, on the two pro- 
totypes, NGC 4395 and POX 52. Interestingly, the SED of 
NGC 4395 differs markedly from those of both quasars and 
typical low-luminosity AGNs (Moran et al. 1999). Specifi- 
cally, the big blue bump, which dominates the SED of luminous 
Seyferts and quasars at optical and UV bands, is absent, lead- 
ing to an optical-to-X-ray spectral index of a ox = -0.97 (De- 
wangan et al. 2008). NGC 4395 is a rare example of a low- 
mass AGN with a relatively low Eddington ratio (Xboi/^Edd) 
of 1.2 x 10~ 3 . Most likely because of strong selection ef- 
fects, optically selected samples of low-mass BHs tend to have 
higher Lboi/^Edd, ranging from < 0.01 to ~ 1 with a median 
of 0.2 (Dong et al. 2012b). POX 52, with L hol = 1.3 x 10 43 
erg s _1 and Lboi/^Edd = 0.2-0.5, exhibits a much more nor- 
mal SED that is broadly similar to that of a scaled-down ver- 
sion of radio-quiet quasars (Thornton et al. 2008). The broad- 
band SED of NGC 405 1 , a narrow-line Seyfert 1 galaxy with 
M BH = (1.7 ±0.5) x 1O 6 M (Denney et al. 2009) accreting at 
1% - 5% of LEdd, is best fit by a relativistically outflowing jet 
model (Maitra et al. 201 1). 

This paper reports a detailed multiwavelength data analysis 
of a low-mass BH hosted in the pseudobulge of UM 625. First 
noted as a "neutral compact spherical disc galaxy" by Zwicky 
et al. (1975), the blue color of UM 625 (B-V = 0.67 mag; 
Salzer et al. 1989) placed it in early catalogs of blue compact 
galaxies (e.g., Campos-Aguilar et al. 1993). A high spatial 
resolution optical image was obtained with HST as part of the 
nearby AGN survey of Malkan et al. (1998), in which UM 625 
was noted as an SO galaxy with a partial nuclear ring structure. 
Later analysis of a near-infrared (NIR) HSTfNICMOS image 



revealed a point-like nucleus embedded in an exponential sur- 
face brightness profile (Quillen et al. 2001; see also Hunt & 
Malkan 2004). An HST ACS/HRC image taken in the near- 
ultraviolet (NUV) shows a very bright, compact, yet partially 
resolved nucleus (Munoz Marin et al. 2007). 

An optical spectrum of UM 625 was first acquired as 
part of the University of Michigan objective-prism survey for 
emission-line galaxies (Mac Alpine & Williams 1981), from 
which its common name originates. Since then UM 625 has 
been classified as a Seyfert 2 galaxy (or occasionally as an H 
II galaxy) in light of its extremely strong [Om] AA4959,5007 
lines and the apparent absence of broad emission lines (e.g., 
Salzer et al. 1989; Terlevich et al. 1991). The spectral classifi- 
cation was re-examined critically by Dessauges-Zavadsky et al. 
(2000); using emission-line measurements from the literature, 
they classified it as a Seyfert 2 according to the three standard 
optical line-ratio diagnostic diagrams of Veilleux & Osterbrock 
(1987). Using a high-quality Sloan Digital Sky Survey (SDSS; 
York et al. 2000) spectrum with a resolution R « 2000, we 
clearly detect both broad Ha and Hf3 and revise the spectral 
classification of UM 625 to a Seyfert 1 . Combining the broad 
line width and luminosity, we estimate a BH mass of 1.6 x 10 6 
M , placing UM 625 in our sample of AGNs with low-mass 
BHs (Dong et al. 2012b). 

We assume a cosmology with Hq = 70 km s _1 Mpc -1 , £l m = 
0.3, and ^ A = 0.7. At a redshift of z = 0.0250, UM 625 has a 
luminosity distance of 109.1 Mpc. 

2. ANALYSIS OF THE OPTICAL SPECTRUM 
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UM 625 was spectroscopically observed by SDSS on 18 June 
2002 UT with an exposure time of 4803 s exposure. It was clas- 
sified as a galaxy by the spectroscopic pipeline of the SDSS 
Fourth Data Release (Adelman-McCarthy et al. 2006). In the 
course of our systematic spectral fitting of all extragalactic ob- 
jects in the SDSS, we first noted UM 625 because it showed 
evident broad Ha and H/3 emission; its inferred virial BH mass 
placed it in the sample of type 1 AGNs with BH masses be- 
low 2 x 10 6 M© described by Dong et al. (2012b). The de- 
tails of the spectral analysis are given in Dong et al. (2012b). 
Here we present a brief description of the continuum modeling 
and emission-line profile fitting, which are based on the MPFIT 
package (Markwardt 2009) that performs x 2 -minimization by 
the Levenberg-Marquardt technique. 

UM 625 has a redshift of z = 0.0250, and its SDSS spectrum 
is dominated by host galaxy starlight. The median signal-to- 
noise (S/N) ratio in the H/3-[Om] and Ha-[Nll]-[Sll] regions 
are 45 and 60 pixel -1 , respectively, high enough to fit accurately 
the continuum and emission lines. We begin by correcting the 
spectrum for Galactic extinction using the extinction map of 
Schlegel et al. (1998) and the reddening curve of Fitzpatrick 
(1999). We model the starlight component with the stellar tem- 
plates of Lu et al. (2006), which were built from the simple 
stellar population spectra (Bruzual & Chariot 2003). The AGN 
continuum is modeled as a power law. The stellar absorption 
lines must be subtracted well to ensure reliable measurement 
of weak emission lines (e.g., Ho et al. 1993, 1997). This is 
achieved by broadening and shifting the starlight templates to 
match the stellar velocity dispersion of the galaxy. As shown 
in Figure 1 (left panel), the fit is good; the absorption features 
are well matched, and the residuals in the emission line-free 
regions are consistent with the noise level. 

Next, we fit the emission lines with Gaussians, using the code 
described in detail in Dong et al. (2005). The spectrum of 
UM 625 is dominated by strong, narrow emission lines such 
as [Om] AA4959,5007, H/3, Ha, and [Nil] AA6548,6583. But 
even a cursory inspection of the spectrum reveals that Ha has 
a clear broad component, even before continuum subtraction; 
after continuum subtraction, a broad component to H/3 also 
emerges. This is reminiscent of the situation of the two pro- 
totypal broad-line AGNs with intermediate-mass BHs, namely 
NGC4395 (Filippenko & Sargent 1989) and POX 52 (Barth et 
al. 2004). Because of their narrowness as well as the high 
S/N ratio of the SDSS spectrum, the doublet lines of both [S II] 
and [Nil] are well isolated. [Sll], [Nil], and the narrow com- 
ponents of the B aimer lines have very similar full widths at 
half-maximum (FWHMs). We fit [Sll] first, assuming that the 
doublet lines have the same width; a good fit is achieved with 
reduced x 2 - 1 • 1 (29 degrees of freedom) by using two Gaus- 
sians for each line. Then we fit the Ha-[Nll] region, taking the 
best-fit model of [S II] as a template to model [N II] and narrow 
Ha. The line ratio of [N II] A6583/A6548 is set to the theoretical 
value 2.96, and their separation is fixed to the laboratory value. 
We use additional Gaussians to model the broad component of 
Ha, starting with one Gaussian and adding in more if the fit 
can be improved significantly according to the F-test. A good 
fit is achieved for broad Ha with just two Gaussians (Figure 1, 
right panel). As the broad H/3 component is weak, we fit the 
total H/3 profile by assuming that its narrow component has the 
same profile as [S II] and that its broad component has the same 
profile as broad Ha. With the exception of [O ill] A4959, 5007, 



Table 1. Fitted emission-line parameters 





Centroid a 


FWHM b 


Flux 


Line 


(A) 


(km s" 1 ) 


(10" 17 erg s" 1 cm" 2 ) 


[Nelll] A3869 


3869.61 ±0.03 


265 ±5 


696 ±11 


[OIII] A4363 


4364.6 ±0.1 


323 ±18 


187 ±8 


Hell A4686 


4687.23 ±0.09 


300 ± 13 


232 ±8 


H/3 (narrow) 


4862.90 ±0.03 


160 


1381 ± 13 


H/3(broad) d 


4864.4 ±2.7 


1801 


317 ±23 


[OIII] A5007 


5008.38 ±0.02 


188 ±4 


7596 ± 39 


[OI] A6300 


6302.36 ±0.04 


197 ±4 


416 ±7 


Hen (narrow) 


6565.09 ±0.04 


160 


5370 ± 38 


Ha (broad) 


6568.0 ±1.3 


1801 ±48 


1660 ± 29 


[Nil] A6583 c 


6585.60 ±0.02 


160 


2291 ±17 


[SII] A6716 


6718.60 ±0.02 


160 ±2 


759 ±9 


[SII] A6731 c 


6732.98 


160 


692 ±9 



a Vacuum rest-frame wavelengths. 
b Corrected for instrumental broadening. 
c Adopting the profile of [S II] A6716. 
d Adopting the profile of broad Ha. 

we fit the other narrow emission lines simply with a single 
Gaussian. For [O III], we modeled each of its doublet lines with 
two Gaussians, one accounting for the bulk component (line 
core) and the other for a weak, yet apparent, blue wing. All the 
line parameters are listed in Table 1 . We make available online 
the data and the detailed fitting parameters 

The ratios of the prominent narrow lines, [O III] A5007/H/3 = 
5.5 and [N II] A6583/Ha = 0.4, place UM 625 within the regime 
of AGNs (Baldwin et al. 1981). With a total H/3 to [O ill] A5007 
ratio of 0.2, UM 625, like NGC 4395 and POX 52, would be 
considered a Seyfert 1.8 in the classification scheme of Os- 
terbrock (1981). From the best-fit model of broad Ha yields 
FWHM = 1801 ±68 krns" 1 and a line dispersion (<Jii ne , the 
second moment of the line profile) of 1583 ± 59 km s _1 (both 
corrected for the SDSS instrumental resolution of 139 km s _1 
FWHM). The luminosity of broad Ha is 2.3 x 10 40 erg s" 1 . 
The observed (not corrected for internal extinction) luminos- 
ity of the narrow Ha component is 3 times higher (7.6 x 10 40 
erg s _1 ), and the luminosity of [O III] A5007 is 1 . 1 x 10 41 erg s _1 . 

To test the reliability of the broad components of B aimer 
lines, we refit the Ha-[Nn] complex and H/3 using as a tem- 
plate for each line the [O III] A5007 profile derived from its 
double-Gaussian model. The results are totally unacceptable; 
the reduced x 2 = 12 and 37, respectively, for the H/3 and the 
Ha-[Nll] complex, with obvious large residuals. This is evi- 
dent by directly comparing the profiles of H/3 and [O III] A5007, 
as shown in the inset of Figure 1 (with [O III] A5007 scaled to 
have the same peak flux density as H/3); H/3 has an additional 
much broader, albeit low-contrast (f\ < 8 x 10~ 17 erg s _1 cm -2 
A -1 ), component that is not present in [Om] A5007. Also evi- 
dent in the inset is that [O III] is slightly broader than H/3 [by ~ 
120 km s" 1 , in height range f\ w (10-20) x 10" 17 erg s" 1 cm" 2 
A -1 ]. We also tried another fitting scheme in which we model 
the narrow B aimer lines with the total profile of [O III] A5007 
and the broad component of the B aimer lines with two Gaus- 
sians; the fitting results are also unacceptable, with reduced 
X 2 = 7 and 23, respectively, for the H/3 and Ha-[Nll] regions. 

The extinction of the broad-line region (BLR) can be derived 



4 Available at http://staff.ustc.edu.cn/~xbdong/Data_Release/IMBH_DR4/ , together with auxiliary code to explain and demonstrate the fitting procedure. 
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Fig. 2.— The archival HST images of UM 625 and 2-D imaging decomposition by GALFIT. From top to bottom are F606W (roughly R band; WFPC2/PC1), 
F330W (U band; ACS/HRC), F160W (H band; NICMOS/NIC1) and composite images of the three bands, respectively. The left column shows the original image, 
the middle column the GALFIT model (PSF + Sersic + disk), and the right column the residual image. All images are oriented with north up and east to the left; the 
red line marks a scale of 1" (~ 0.5 kpc). The green polygon denotes the star-forming ring region, which has been masked out in the GLAFIT fitting. 



from the observed Balmer decrement Ha/H/3 for normal AGNs 
(Dong et al. 2008). Assuming the extinction curve of the Small 
Magellanic Cloud (Hopkins et al. 2004; Wang et al. 2005) and 
an intrinsic broad-line Ha/H/3 = 3.1 (Dong et al. 2008), we 
get E(B-V) = 0.57 mag for the BLR. Likewise, the observed 
narrow Ha and H/3 gives E(B-V) = 0.25 mag for the NLR. 

3. ANALYSIS OF THE IMAGES 

There are archival HST images in three bands for UM 625: 
F330W (roughly the U band of the Johnson system), F606W 



(/?), and F160W (#), observed through ACS, WFPC2, and 
NICMOS, respectively. From these images, we can see clearly 
that there is a bright point-like source present in the center of 
UM 625, and that the galaxy is almost round in shape and has 
no spiral arms, indicating an elliptical/spheroidal or a face-on 
SO galaxy. Closer inspection reveals that there is, in addition, 
a nuclear semi-ring on scales of ~ 150-400 pc; this feature, 
most conspicuous in the F330W image, is resolved into sev- 
eral knots, presumably sites of intense star formation. The ring 
region is indicated in Figure 2 with the green polygon. These 



Jiang et al. 




18 


. i i i i | i i i i | i i i i | i i i i . 




■ \ 


_ 20 


- \ 


U 


\ 


CD 


• # 


en 


• 


e 22 


^""\ ""• 


D 


^w*. 


cn 


^^^^ 


£ 24 


: ^^^v, : 


i 


\, 


26 


- ^^>->L 


28 





10 15 

smo (arcsec) 



20 



Fig. 3. — (Left) The SDSS three-color (gri) composite image of UM 625; the yellow line marks a scale of 10". (Right) The r-band surface brightness profile 
(dzlcr error bars) as a function of semi-major axis radius, with sky subtracted. The red dotted line denotes the limiting surface brightness of 26.57 mag arcsec -2 . We 
fit the profile at r > 10" with an exponential disk model (blue line), which gives a scale length of 4"21 (2.12 kpc). 



high-resolution images give us a distinct opportunity to study 
not only the multiwavelength properties of the AGN, but also 
its host galaxy. The HST images, however, suffer from one ma- 
jor limitation: their field of view is too small to properly sample 
the disk component or to reliably measure the sky background. 
Fortunately, UM 625 was imaged with SDSS. Thus, before 
we discuss the HST images in detail, we first turn to the SDSS 
images. Our strategy is to use the SDSS r-band image, which 
closely approximates the HST F606W band, to constrain the 
photometric parameters of the disk component, which is other- 
wise difficult to determine from the HST data alone. With the 
disk thus constrained, we will use the HST F606W image to 
measure the parameters of the AGN point source and the bulge, 
which are central to our scientific analysis. The F330W and 
F160W images provide further photometric points for the nu- 
clear SED, as well as color information to diagnose the stellar 
population of the bulge. 

3.1. SDSS Imaging 

UM 625 was observed by SDSS in ugriz on 24 May 2001 
UT (Figure 3). Although the standard SDSS images have a rel- 
atively short exposure time of only 54 s per filter, the drift- scan 
mode in which the imaging survey was conducted (Gunn et al. 
1998) ensures very accurate flat-fielding. Moreover, for most 
galaxies the field of view is quite large. Both factors are crucial 
for obtaining an accurate measurement of the sky background 
and its associated error, which determines the limiting surface 
brightness sensitivity. For most galaxies the azimuthally aver- 
aged surface brightness profile can be measured reliably down 
to n r « 27 mag arcsec -2 , which is deep enough to study the 
outer structure of galaxies (Pohlen & Trujillo 2006; Erwin et al. 
2008; Li etal. 2011). 

We chose to work with the r band, which is closest to the 
HST F606W filter. We first mask all photometric objects in the 
field identified either by Sextractor (Bertin & Arnouts 1996) 

5 IRAF (Image Reduction and Analysis Facility) is distributed by the National 
agreement with the National Science Foundation. 

^http: //www, astro . yale ."e du/dokkum/lacos mic/| 



or by the SDSS photometric pipeline and then run the IRAI0 
task ellipse to fit isophotes with a linear step of 2 pixels 
between successive steps, allowing the center, position angle, 
and ellipticity of each ellipse to vary. The sky background level 
is determined from the best-fit ellipses that have constant sur- 
face brightness with respect to radius, and the uncertainty on 
the background is estimated from the root mean square fluctua- 
tions (a) about the mean value. The limiting surface brightness, 
defined as 3 a, is 26.57 mag arcsec -2 . After sky subtraction, the 
azimuthally averaged one-dimensional (1-D) surface brightness 
profile is extracted in logarithmic steps, to increase the S/N in 
the noisier outer regions. We can clearly see that the profile 
(right panel of Figure 3) at radii larger than ~ 10" behaves like 
an exponential disk, the scale length of which is 4. // 21±0. // 17 
(2.12 ± 0.09 kpc). We do not decompose the inner regions of 
the galaxy with this data set, as the SDSS image lacks sufficient 
resolution; for that, we turn to the high-resolution HST images 
below. 

3.2. HST: Optical Image 

The F606W image is our best choice for probing the host 
galaxy structure because (1) the bandpass is relatively red and 
thus less sensitive to dust and young stars, (2) the field of view 
of WFPC2/PC1 is larger than that of ACS/HRC and NICMOS, 
and (3) its point-spread function (PSF) has less extended wings 
that NICMOS. A single 500 s exposure was obtained on 21 July 
1994 (Proposal ID: 5479), with UM 625 placed near the cen- 
ter of the Planetary Camera detector (PCI), which has a plate 
scale of 0."046 pixel -1 . We remove cosmic rays using LA Cos- 
miqj (van Dokkum 2001), which can detect cosmic ray hits of 
arbitrary shape and size. Two of the pixels in the center of the 
galaxy are saturated; we masked them out in the analysis below. 

Precise sky subtraction is of great importance to obtain accu- 
rate photometric measurements and structural decomposition. 
From the surface brightness profile of the SDSS r-band image, 

Optical Astronomy Observatory, which is operated by AURA, Inc., under cooperative 
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Fig. 4. — Radial surface brightness profile of UM 625 extracted from 
the HST WFPC2/PC1 F606W image (red crosses). The profile with the star- 
forming ring masked is displayed with black dots (with ±1<j error bars). For 
comparison we also plot the profile of the PSF created by Tiny Tim (blue dotted 
line) and the PSF profile from a star taken from the WFPC2 PSF library (or- 
ange diamonds). The PSF profiles are scaled to have the same central surface 
brightness as UM 625. The green squares represent the SDSS r-band profile 
with sky subtracted. 

we know that UM 625 extends close to a radius of 20", which 
means that it fills nearly the entire field of PCI (radius ~ 18"). 
We estimate the background level and its uncertainty from the 
outermost edges of the PCI chip and from the flanking WF2 
chips; the limiting surface brightness of the F606W image is 
25.31 mag arcsec -2 . To examine the effect of the star- forming 
ring on the underlying galaxy structure, we also plot the surface 
brightness profile with the ring region (the region enclosed by 
the green polygon as indicated in Figure 2) masked out. Fig- 
ure 4 shows that the circumnuclear star- forming region mainly 
affects the profile on scales of 0."2< r <0."8; the outer profile 
(r > 2") agrees well with the SDSS r-band profile. 

Separating the central AGN light from the host galaxy 
starlight requires knowing the PSF to a high accuracy. Unfor- 
tunately there is no bright star in the field. According to the 
image simulations of Kim et al. (2008), a synthetic PSF gen- 
erated by the Tiny Tim software (Krist 1995) works reasonably 
well. We have also searched for empirical stellar PSFs from the 
WFPC2 PSF Librar)Q and chosen the one closest to the loca- 
tion of UM 625 on PCI. As shown in Figure 4, the stellar PSF 
agrees well with the Tiny Tim PSF, except in the outer regions 
where the empirical PSF is much noisier. For the rest of the 
analysis, we simply adopt the Tiny Tim PSF. 

After background subtraction, we perform a two- 
dimensional (2-D) decomposition of UM 625 using GALFIT 
(Peng et al. 2002, 2010). The AGN is represented by a point 
source modeled with the Tiny Tim PSF, and the galaxy is mod- 
eled by bulge, fit with a Sersic (1968) r l l n function, and a 
disk, fit with an exponential function (equivalent to n = 1). A 
single-component model for the host leaves unacceptably large 

7 http://www.stsci.edu/hst/wfpc2/software/wfpc2-psf-form.html 
8 http ://drizzlepac . stsci.edu/ 
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Fig. 5. — One-dimensional representation of the three-component 2-D 
GALFIT model applied to the WFPC2/F606W image of UM 625: PSF for the 
nucleus (black dotted line), n = 1.60 Sersic function for the bulge (red dashed 
line), and an exponential function for the disk (blue dot-dashed line). The sum 
of the three components is shown as the black solid line. The observed data 
are plotted as black symbols with ±1<t error bars. The bottom panel shows the 
residuals between the data and the best-fit model. 

residuals. The best-fitting two-component model (Figure 5) 
gives a central component with Sersic index n- 1.60 and ef- 
fective radius r e = l."37 (693 pc), and a disk with scale length 
of 3. // 86 (1.95 kpc), consistent with the results from the 1-D de- 
composition of the SDSS r-band image. The bulge-like compo- 
nent is mildly disky, with c = -0. 1 1 . The best-fit parameters are 
summarized in Table 2. We identify the central n = 1 .60 compo- 
nent with a pseudobulge, because pseudobulges generally have 
n < 2 (Kormendy & Kennicutt 2004). The residual image (top 
right panel of Figure 2) reveals a faint ring-shaped structure on 
the eastern side of the galaxy, opposite to the masked ring re- 
gion. 



3.3. HST: NUV Image 

The F330W image was observed on 21 March 2003 using 
ACS/HRC as part of a study on the starburst-AGN connec- 
tion (Proposal ID: 9379). ACS is located away from the op- 
tical axis of HST, and so it suffers from significant geometric 
distortion that is not corrected by the internal optics. For ease 
of rejecting cosmic rays, the total exposure was divided into 
two equal exposures of 10 min each. These images are then 
combined using astrodrizzleQ, a new software replacing 
mult idr i z z le, to remove cosmic ray hits and to correct for 
the geometric distortion. The pixel scale of the combined image 
is 0. // 025. A narrow ring-like star- forming region, ornamented 
with several bright knots, encircles nearly half the nucleus. We 
begin with creating two synthetic PSFs for each exposure using 
Tiny Tim. They are centered in the same position as the nucleus 



Table 2 . GALFIT fiffiJ&gtSlion. 



Band 


Component 
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M 


n 


r C'/pc) 


b/a 


c 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


F606W 


PSF 


17.89 


-17.54 












Sersic 


15.97 


-19.62 


1.60 


1.37/693 


0.87 


-0.11 




Exp Disk 


16.67 


-18.74 


[1.0] 


3.86/1948 


0.79 






Ring 


18.20|18.16 


— 17.21| — 17.25 










F330W 


PSF 


17.88 


-17.14 












Sersic 


16.98 


-18.30 


2.64 


1.02/515 


1.00 


-0.55 




Exp Disk 


17.73 


-17.66 


[1.0] 


[3.86/1948] 


0.55 






Ring 


17.64|17.62 


-17.52|-17.54 










F160W 


PSF 


16.77 


-18.49 












Sersic 


13.68 


-21.54 


1.54 


0.98/491 


0.82 


-0.06 




Exp Disk 


[14.56] 


[-20.74] 


[1.0] 


[3.86/1948] 


[0.79] 






Ring 


17.43|17.47 


-17.79|-17.75 











Note. — Col. (1): HST filter. Col. (2): Components used in the fitting schemes. Col. 
(3): The integrated magnitudes on the Vega system, not corrected for Galactic extinction. 
The ring magnitude is derived from the residual images, with two values given corresponding 
to the magnitude in the ring region and aperture photometry, respectively. Col. (4): The 
absolute Johnson R, U and H magnitude after Galactic extinction correction. We assume 
a continuum (f u oc z^ -0-3 , which is the same as our optical spectral fitting) for the central 
AGN, an Sb galaxy for the Sersic component, and an Sc galaxy for the disk and ring, using 
templates from Kinney et al (1996). The SWIRE Template Library (Polletta et al. 2007) 
is used to convert F160W to H magnitude. Col. (5): The Sersic index. Col. (6): The 
effective radius of the Sersic component or scale length of exponential disk, in units of arcsec 
and pc, respectively. Col. (7): Axis ratio. Col. (8): Diskiness (negative) /boxiness (positive) 
parameter, defined in Eqn. (3) of Peng et al. (2002). The brackets mean that they are fixed. 
The formal errors given by GALFIT are all tiny: < 0.05 for magnitude and Sersic index, 
< 0."1 for r. 



in the HRC image to properly reproduce the geometric changes 
involved in the processing. During the modeling, we have also 
considered the dofocus offset (~ 6 fim) given by web-based 
modefl The offset is mostly due to spacecraft breathing ef- 
fects. The two Tiny Tim PSFs, still distorted, are almost iden- 
tical except for slightly differences due to defocus offsets. We 
combined them with astrodrizzleinthe same manner as 
the science images, producing a distortion-corrected PSF im- 
age. 

We decompose the image using GALFIT with the region 
containing the star-forming ring masked out. Due to the small 
field of view of the ACS/HRC (26" x 29"), we cannot indepen- 
dently measure the sky background; we set it as a free param- 
eter in the fit. Nor can be obtain any meaningful, independent 
constraint on the parameters of the disk. Our initial attempt to 
fit the image with a PSF + Sersic + disk model failed to con- 
verge. As our primary goal for the F330W image is to mea- 
sure the brightness of the nucleus, we simply set the disk scale 
length to be identical to that obtained from the F606W fit. The 
resulting fit yields a bulge with n = 2.64 and r e = 1 "02 (515 pc), 
and a point source flux that is very insensitive to the assump- 
tions of the disk component. The AGN magnitude changes by 
only ±0.01 mag depending on whether the disk is included or 
not. If we further force the bulge to have the same n and r e as 
derived from the F606W fit (i.e. allowing only the luminosity 
to adjust), the resulting AGN magnitude changes by 0.06 mag. 



the standard pipeline calnica within IRAF/STSDAS. This 
task corrects for the nonlinearity of the detector and removes 
the bias value, dark current, amplifier glow, and shading; how- 
ever, it does not remove the ghost "pedestal" effect produced 
by the variable quadrant bias (see, e.g., Hunt & Malkan, 2004). 
To correct for this effect, we apply the pedsub task to the 
calibrated image, quadrant by quadrant, determine the shifts 
of the dithered images with xregister, register them with 
imshif t, and finally combine them with imcombine. 

UM 625 extends far beyond the ~ 1 1" x 1 1" field of view of 
NIC1. In spite of this, a GALFIT model consisting of a PSF 
+ Sersic bulge can still give a fairly robust measurement of the 
nuclear magnitude. During the fitting, we set the background 
free and masked the ring region as before. It proved to be im- 
possible to independently measure the disk component at all. 
Since we know from the optical images that a disk is present, we 
constrain it by fixing its structural parameters to those derived 
from the F606W image and normalizing its luminosity such that 
R-H = 2 mag, a value typical of disks (e.g., MacArthur et al. 
2004). However, we note that none of the parameters for the 
bulge or the AGN point source are significantly affected by our 
assumptions for the disk component. The best-fit model 
yields n = 1.54 and r e = 0."98 (491 pc) for the bulge. If we fix 
the bulge n and r e to the values derived from the F606W fit, the 
AGN magnitude changes by 0.21 mag. 

4. OTHER MULTIWAVELENGTH DATA 



4.1. X-ray Observations 
UM 625 was observed by XMM-Newton on three sepa- 



3.4. HST: NIR Image 

The NICMOS/NIC1 F160W image was observed on 31 July 
1997 (Proposal ID: 7328), in three 256 s exposures dithered rate occasions between 2004 and 2008 (ObsID: 0200430901, 

in an "L"-shaped pattern. Each exposure was processed with 0505930101, 0505930401) and once by Chandra in 2008 (Ob- 

9 http://www.stsci.edu/hst/observatory/focus/FocusModel. We did not include the offset information for the F606W image because it is only available for data 
taken after 2003. 
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Fig. 6. — X-ray spectra of UM 625. The upper panel shows the data and best-fit model; the bottom panel shows the ratio between the data and model. The dotted 
lines represent the two power-law components used in the fitting convolved with the response of the detectors. Four data sets are shown: Chandra data from 2008 
(black) and XMM-Newton PN data from 2004 (blue), 2007 (red), and 2008 (green). 



sID: 9557). We used XMM-Newton SAS version 12.0.1 and the 
calibration file of June 2012 to reduce the EPIC PN (Struder et 
al. 2001) archival data. After filtering background flares, the 
net exposure times are 6.4 ks, 13 ks, and 13 ks, respectively, for 
the three PN observations. The source was detected at a large 
offset angle (~ 10" -11") from the field center for the later two 
PN observations. To extract the source and background spectra, 
we defined the source region as a circle with radius ~ 30" -40" 
centered at the source position and the background regions as 
three or four circles, far removed from the CCD edges and other 
sources, with radii ~ 30" - 40" around the source region. The 
Chandra ACIS-S (Garmire et al. 2003) data were reduced with 
CIAO 4.4.1 (Fruscione et al. 2006) and the CALDB 4.4.10 
database, adopting standard procedures. The net exposure time 
was 49 ks. Because of the large off-axis angle (10."6), we ex- 
tracted the spectra using a 20"-radius circle for the source and 
four circles around the source region for the background. All 
XMM-Newton and Chandra spectra were rebinned to have a 
minimum of 5 counts in each energy bin after background sub- 
traction, and they were fit simultaneously with the same models 
using Xspec (Arnaud 1996). In view of the low source counts, 
we fit the spectra using the Cash- statistic instead of x 2 - 

We first fit the spectra over the energy range 0.2-10 keV 
using a single absorbed power-law model with the minimum 
absorption set to the Galactic value of Nu = 4.04 x 10 20 cm~ 2 
(Kalberla et al. 2005). The fit yields a column density equal 



to the Galactic value, a photon index Ti = 2.65^o;i5 ? and C- 
statistic = 145 (125 degrees of freedom). The residuals, how- 
ever, clearly show a flat, hard excess above 3 keV, which im- 
plies that an obscured component exists. We then added an 
absorbed power-law component with a photon index fixed to 
T 2 = 1.90, typical of AGNs (Shemmer et al. 2008; Constantin 
et al. 2009). As shown in Figure 6, the fit is improved signifi- 
cantly, resulting in C-statistic = 116 (123 degrees of freedom). 
The soft-band power law, with T\ = 2.86^;}^, is slightly steeper 
than the previous fit. The absorption column density of the ab- 
sorbed power-law component is Nu = 9.7^ 8 8 x 10 22 cm~ 2 . The 
observed 0.5-2 keV and 2-10 keV fluxes are 2.5!°; 2 x 10~ 14 and 
4.6+Q4 x 10" 14 erg s" 1 cm -2 , corresponding to 4.0+[j; 4 x 10 40 and 



-0.4 

6.5^6 x !0 4 ° erg s~ l 



respectively. Assuming that the obscured 
power-law component is the intrinsic emission from the corona, 
the unabsorbed 2-10 keV flux is 6.7 x 10~ 14 erg s _1 cm -2 , cor- 
responding to 9.5 x 10 40 erg s _1 . The monochromatic flux at 
2 keV is 7.8 x 10~ 32 erg s _1 cm -2 Hz -1 . As the spectra are cut 
off above 10 keV and no variations are detected in each single 
exposure or among the four exposures, we cannot exclude the 
possibility that the absorption is Compton-thick. 

If we link the spectral slope of the hard component to the soft 
one, we obtain T = 2.84 ± 0.16, N H = 14.6! 23 8 x 10 42 cnT 2 , and 
a scattering fraction of 10% if we attribute the soft component 
to scattered nuclear emission. We note that this fraction is con- 
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si stent with those commonly reported in other obscured AGNs 
(0.1%- 10%; Noguchi et al . 2010). 

4.2. Ultraviolet Observations 

The Optical Monitor (OM) on XMM-Newton took data for 
UM 625 in the UVW1 (2910 A) and UVW2 (2120 A) fil- 
ters, each for 1 ks, during the 2004 X-ray observations. Be- 
cause of the coarse resolution of the OM (PSF « l."8), the 
nucleus is not resolved from the host galaxy, and we consider 
these global (AGN plus host) flux measurements. We adopt 
an aperture radius of 10" to compute the source flux after es- 
timating the sky background from the region between a radius 
of 15" and 20". The flux is corrected for Galactic extinction 
using the maps of Schlegel et al. (1998) and the reddening 
curve of Fitzpatrick (1999), resulting in ^(2910 A)=(4.5 ± 
0.3) x 10~ 27 erg s~ l cm" 2 Hz" 1 and ^(2 120 A)=(2.7 ± 0.2) x 
10~ 27 erg s _1 cm -2 Hz -1 . 

Apart from the XMM-Newton OM data, UM 625 was also 
imaged simultaneously in the near-UV (NUV; 2316 A) and far- 
UV (FUV; 1539 A) bands of GALEX during its All-sky Imag- 
ing Survey (AIS), for a total exposure time of 1 12 s on 6 April 
2004. UM 625 was later reobserved during the Medium Imag- 
ing Survey (MIS) on 18 May 2009, for a total of 1578 s; we 
adopt the MIS data (Bianchi et al. 2012). After correcting for 
Galactic extinction, we obtain an NUV flux density of 7^(2316 
A)=(3.17±0.04) x 10~ 27 erg s~ l cm" 2 Hz" 1 and an FUV flux 
density of ^(1539 A)=(1.42±0.04) x 10~ 27 erg s~ l cm" 2 Hz" 1 . 
The NUV flux density is slightly larger than that of the OM 
UVW2 band, whose effective wavelength is only ~ 200 A 
shorter. 

4.3. Radio Observations 

UM 625 was detected by Faint Images of the Radio Sky at 
Twenty-cm (FIRST; Becker et al. 1995) using the Very Large 
Array in its B configuration. The peak and integrated 20 cm 
flux density from the FIRST catalog (White et al. 1997) is 
1.89 and 2.69 mJy, corresponding to a monochromatic radio lu- 
minosity at 20 cm of 2.66 x 10 28 and 3.79 x 10 28 erg s" 1 Hz" 1 . 
These measurements are derived by fitting a 2-D Gaussian func- 
tion to the source, using a map generated from twelve coad- 
ded images adjacent to the pointing center; the map has l."8 
pixel -1 , a resolution of FWHM = 5", and an rms noise of 0.152 
mJy beam -1 . Customarily, the radio-loudness parameter R is 
defined as the ratio of flux densities between 6 cm and 4400 A. 
Assuming a radio spectral index a r = -0.46 (Lai & Ho 2010), 
f y {6 cm) = 1.55 x 10 -26 erg s -1 cm -2 Hz -1 , which, in combi- 
nation with 7^(4400 A) determined from the SDSS spectral fit- 
ting, gives R = 5.0. Although UM 625 is formally radio-quiet 
according to the widely used division of R = 10 for radio-loud 
and radio-quiet AGNs (Kellermann et al. 1989), its level of ra- 
dio activity is still somewhat unusual for low-mass AGNs (e.g., 
Greene et al. 2006). 

4.4. Infrared Observations 

UM 625 is also contained in the Infrared Astronomical Satel- 
lite (IRAS) faint source catalog (Moshir et al. 1990). The faint 
source catalog contains data for point sources in unconfused re- 
gions with flux densities typically above 0.2 Jy at 12, 25, and 
60 /im, and above 1.0 Jy at 100/im. UM 625 is reliably de- 
tected at 60 /im with a flux density of 0.32 ± 0.08 Jy, but only 

10 http://sundog. stsci.edu/cgi-bin/searchfirst 



upper limits are given for the other three bands (0.13, 0.12, and 
0.77 Jy for 12, 25, and 60 /im, respectively). 

Wide-field Infrared Survey Explorer (WISE; Wright et al. 
2010), the mission most comparable to IRAS yet with a sen- 
sitivity more than 100 times higher at 12/im, has mapped the 
entire sky in four bands centered at 3.4, 4.6, 12, and 22 /im. 
UM 625 is detected with high S/N in all four bands, with a flux 
density, converted from the magnitudes in the WISE All- sky 
Source Catalog, of 4.1, 4.2, 24, and 95 mJy, respectively. 

For completeness, we take the following NIR magnitudes 
from the Two Micron All Sky Survey point source catalog 
(Skrutskie et al. 2006): J= 14.25 ±0.05, H= 13.60 ±0.06 and 
Kg = 13.24 ± 0.05 mag, which correspond to a flux density of 
3.3, 3.8, and 3.4 mJy, respectively. 

5. RESULTS AND DISCUSSION 

5.1. Black Hole Mass 

With the detection of broad emission lines, we can estimate 
the mass of the central BH using commonly used virial mass 
estimators for broad-line AGNs. We calculate the virial BH 
mass using the Ha formalism given in Xiao et al. (2011; their 
Equation 6), which is based on Greene & Ho (2005b, 2007b) 
but updated with the more recent relationship between BLR 
size and luminosity of Bentz et al. (2009). For FWHM(Ha) 
= 1801 km s _1 and a broad Ha luminosity of 2.4 x 10 40 erg s" 1 , 
Mbh = 1.6 x 10 6 Mq, which justifies inclusion in the low-mass 
BH sample of Dong et al. (2012b). 

The uncertainty of the above BH mass estimate is not well 
understood. As pointed out by Vestergaard & Peterson (2006), 
the statistical accuracy of the masses from single-epoch virial 
mass estimators is a factor of ~ 4 (la), and, for individual mass 
estimates, the uncertainty can be as large as an order of mag- 
nitude. We note that recently Wang et al. (2009) recalibrated 
the BH mass formulas based on single-epoch spectra, stressing 
the nonlinear relation between the virial velocity of the BLR 
clouds and the FWHM of single-epoch broad emission lines. 
This nonlinearity probably arises from several kinds of non- 
virial components incorporated into the total profile of broad 
emission lines in single-epoch spectra (see §4.2 of Wang et al. 
[2009] for a detailed discussion, as well as Collin et al. [2006] 
and Sulentic et al. [2006]). The formulas of Wang et al. (2009), 
calibrated using reverberation mapping data available to date, 
which span the mass range from M B h ~ 10 7 to 10 9 M , does 
not cover the low-mass regime of interest in this paper. Further- 
more, we would like to point out that nonvirial components, if 
any, may be less significant in UM 625 than in other AGNs be- 
cause its broad Ha profile is roughly symmetrical and close to 
a Gaussian. 

5.2. Host Galaxy 

The abundant archival HST images covering a wide wave- 
length baseline provide us an excellent opportunity to study the 
host galaxy of UM 625 with the central AGN point source re- 
moved. Our 2-D decomposition of the optical (F606W, R band) 
WFPC2/PC1 shows that the stellar distribution consists of two 
main components: (1) a dominant compact, bulge-like compo- 
nent with Sersic n = 1.60, r e = 693 pc, and M R = -19.62 mag; 
(2) an extended disk (n « 1) component with a scale length of 
~ 2 kpc. The disk is more evident in the SDSS r-band im- 
age owing to its deeper limiting surface brightness and larger 
field of view, but its photometric parameters are well recovered 
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with the WFPC2/PC1 image alone. The very limited field of 
view of the ACS/HRC and NICMOS/NIC1 images prevent us 
from placing any meaningful constraints on the disk in the UV 
and NIR. We suspect that even the structural parameters of the 
bulge may be partly compromised in these bands as a result of 
the small field size. 

Taking the F606W decomposition as reference, the bulge-to- 
total light ratio (B/T) in the R band is 0.66. If we use a total 
host luminosity calculated from the SDSS r-band image (after 
subtracting the AGN point source derived from the H ST decom- 
position), which may be more reliable because of its larger field 
of view, the bulge-to-total ratio reduces slightly to B/T = 0.60. 
Either of these values lies within the range of B/T for SO galax- 
ies (e.g., Simien & de Vaucouleurs 1986), in agreement with 
the absence of spiral arms in the main disk of the galaxy. The 
star-forming ring, which resembles a tightly wound spiral, in 
UM 625 is located in the circumnuclear region (~ 0.5 kpc) re- 
solved only by HST observations. Interestingly, the UV-optical 
and optical-NIR colors of the bulge of UM 625 are much bluer 
than those of typical SO galaxies; they more closely resemble 
those of an Sb spiral. This is apparent from comparing the ob- 
served F330W-F606W and F606W-F160W colors with syn- 
thetic colors calculated using calcphot in IRAF/SYNPHOT 
package for various galaxy templates (Kinney et al. 1996; Pol- 
letta et al. 2007). 

In view of its blue color, low Sersic index (n < 2), and disky 
isophote shape (c < 0), the bulge of UM 625 can be convinc- 
ingly categorized as a pseudobulge, commonly present in low- 
luminosity disk galaxies, including some SOs (see Kormendy & 
Kennicutt 2004 for a review). The presence of an ongoing cen- 
tral star formation, manifested through the nuclear star- forming 
ring (Section 5.3), further supports the pseudobulge interpreta- 
tion. This is consistent with previous studies, which find that 
the majority of the host galaxies of low-mass BHs with disks 
are likely to contain pseudobulges rather than classical bulges 
(Greene et al. 2008; Jiang et al. 2011b). The only possible 
anomaly is that UM 625 has a much higher B/T than usual; for 
example, the objects in the sample of Jiang et al. (201 lb) with 
a detected disk component have an average B/T = 0.23. 

5.3. Nuclear Star-forming Ring 

As mentioned above, the circumnuclear region of UM 625 
contains a blue (F330W-F606W = -0.73 mag, corrected for 
Galactic extinction) semi-ring on scales of ~ 150-400 pc. It 
is most prominent in the NUV (F330W) band, which shows a 
number of bright, compact knots, reminiscent of star- forming 
galaxies with nuclear hotspots (e.g., Barth et al. 1995). As- 
suming that the UV light mainly arises from young stars, we 
can use the integrated luminosity to estimate the star forma- 
tion rate. Summing the flux of the residuals above the smooth 
galaxy model within the ring region (green polygon in Figure 2) 
yields a luminosity of L v = 2.0 x 10 27 erg s _1 Hz -1 after correct- 
ing for Galactic extinction, which corresponds to a star forma- 
tion rate of 0.28 M yr" 1 (Kennicutt 1998; his Equation 1). As a 
check, aperture photometry between radius 0."15 and \" in the 
residual image yields nearly the same flux, only 0.02 magnitude 
higher. 

Given that the central region of UM 625 experiences signifi- 
cant ongoing star formation, as evidenced by the UV light, it is 
interesting to note that some of the narrow emission lines, es- 
pecially the B aimer lines, are inevitably contaminated by stellar 
photoionization. From Equation (2) of Kennicutt (1998), a star 



formation rate of ~ 0.3 M yr -1 produces an Ha luminosity of 

3.5 x 10 40 erg s _1 , which is roughly half of the total observed 
narrow Ha emission. This is nonnegligible. 

The star-forming region can potentially contaminate the 
emission in the radio and X-ray bands as well, due to the con- 
tribution from high-mass X-ray binaries, young supernova rem- 
nants, and hot interstellar plasma. According to the empirically 
calibrated linear relation between star formation rate and X-ray 
luminosity of Ranalli et al. (2003), the ring region produces 
£0.5-2 kev = 1.3 x 10 39 erg s _1 ; this is only a few percent of 
the observed X-ray luminosity, and so the AGN totally dom- 
inates the X-ray emission. By contrast, the Ranalli et al.'s re- 
lation between star formation rate and radio emission predicts 
L\a GHz = 1.1 x 10 28 erg s _1 Hz -1 , which is nearly 30% of the 
integrated monochromatic radio power detected by FIRST. 

5.4. X-ray Spectral Properties 

Although the detection of broad Ha and H/3 qualifies UM 
625 as a type 1 AGN, its X-ray spectrum indicates that it con- 
tains a significant intrinsic absorbing column density of Nu = 
9.7^ 8 8 x 10 22 cm~ 2 . This apparent disagreement between opti- 
cal and X-ray classification is not uncommon (e.g., Garcet et al. 
2007), as it can arise if the gas responsible for the X-ray absorp- 
tion is highly ionized, instead of neutral, so that the accompa- 
nying dust would sublimate to yield a much smaller dust-to-gas 
ratio. Indeed, such absorbed type 1 AGNs are usually charac- 
terized by either complex or warm/ionized absorption (Malizia 
et al. 2012) arising from ionized gas possibly associated with 
a disk wind (Murray et al. 1995) or ionization cones as seen in 
some objects. 

The ratio T between 2-10 keV luminosity to extinction- 
corrected [O III] luminosity is also a powerful diagnostic of nu- 
clear X-ray obscuration. Previous studies have suggested that 
objects with T ^ 0.1 are invariably Compton-thick, whereas 
objects with T ^ 1 are almost exclusively Compton-thin or un- 
obscured (Guainazzi et al. 2005). With observed L2-iokeV ~ 

4.6 x 10 40 erg s _1 and an extinction-corrected L\o nrj = 2.3 x 
10 41 erg s _1 , T - 0.2, placing UM 625 intermediate between 
Compton-thick and Compton-thin. 

The power-law component for the soft X-ray band is likely 
different from the soft X-ray excess commonly seen in type 1 
AGNs since it should have been largely obscured, if present. 
It could arise from a scattered component from the nucleus or 
contamination from the host galaxy; the contribution from the 
star- forming ring is negligible (Section 5.3). For AGNs with 
low-mass BHs, thermal emission from the accretion disk can 
also contribute significantly to the soft X-ray band (e.g., Thorn- 
ton et al. 2008; Miniutti et al. 2009). We have tried to fit our 
X-ray data using a (disk) blackbody model, but the data qual- 
ity is insufficient to reach meaningful conclusions regarding the 
nature of the soft X-rays. 

5.5. Spectral Energy Distribution and Bolometric Luminosity 

We combine all of the photometric data of UM 625 described 
in Sections 3 and 4 to construct its broad-band SED (Figure 7); 
apart from NGC 4395 and POX 52, this is one of the most com- 
plete SEDs available for low-mass AGNs. The median SEDs of 
radio-quiet and radio-loud quasars (Elvis et al. 1994), scaled to 
UM 625 in the optical band, are overplotted for comparison. 

Even after correcting for absorption, the X-ray emission 
of UM 625 is still somewhat weak compared with the 
median SED of quasars. The optical-to-X-ray slope a ox 
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Fig. 7. — The SED of UM 625. The nuclear measurements pertaining to the AGN are plotted in black dots: radio data from FIRST (with the contribution from 
star-forming ring subtracted), NUV, optical, and NIR data from HST decomposition (error bars denote the range of magnitudes allowed by the different GALFIT 
fitting schemes), and unabsorbed X-ray spectrum from Chandra and XMM-Newton (la uncertainty range given by the dot-dashed line). Integrated measurements 
for the entire host galaxy (plus AGN) are plotted with various red symbols. We overplot the median SED of radio-loud (black dotted line) and radio-quiet (black 
dashed line) quasars (Elvis et al. 1994), scaled to the HST nuclear optical point. We also plot the SED of a typical Sb galaxy (green line) from the SWIRE template 
library of Polletta et al. (2007), scaled to the HST optical point for the host galaxy. 



is -1.72, where we adopt the standard definition a ox = 
-03841og[/ I X2500A)// I/ (2keV)] (Tananbaum et al. 1979) and 
estimate ^(2500 A) from ^(5100 A) assuming an optical-UV 
continuum spectral index of -0.44 (Vanden Berk et al. 200 lfl 
By contrast, previous X-ray studies of low-mass AGNs indi- 
cate that a ox is, on average, larger than in high-mass AGNs 
(e.g., Greene & Ho 2007a; Miniutti et al. 2009; Desroches 
et al. 2009; Dong et al. 2012a), falling systematically below 
the low-luminosity extension of the a O x-^(2500A) relation of 
Steffen et al. (2006). The previous samples span a wide range 
in a ox , from « -1.7 to -1. Thus, UM 625 lies among the weak- 
est X-ray sources with the lowest a ox . The origin of the X-ray 
weakness is not known. Some may be simply highly absorbed, 
but others may be intrinsically X-ray weak (Dong et al. 2012a). 
UM 625 may belong to the latter category. 

Low-mass AGNs appear to be exceptionally radio-quiet, if 
not radio- silent. In a Very Large Array 6 cm survey of the 
19 low-mass AGNs from Greene & Ho (2004), Greene et al. 
(2006) detected radio emission from only one (~ 5%), which 
has 7?=2.8. This detection rate is approximately the same as 
that found in the larger sample of Greene & Ho (2007b), whose 
detected sources have R « 1-80. UM 625 is clearly detected 



at 20 cm, at R = 5.0. Even after correcting for possible contam- 
ination from star formation (Section 5.3), the source still has 
R w 3.3. While UM 625 is, strictly speaking, not radio-loud, it 
is still somewhat unusual in that it belongs to the minority of 
low-mass AGNs that show any radio emission at all. 

Despite the fact that the SED of UM 625 still contains many 
gaps in wavelength coverage, it should still offer a more reli- 
able measurement of the bolometric luminosity of the system 
than any estimate based on a single band. Integrating the me- 
dian radio-quiet quasar SED of Elvis et al. (1994) after scal- 
ing it to the nuclear optical point derived from HST for UM 
625, we obtain Lb i = 2.4 x 10 43 erg s _1 , which corresponds to 
^boi/^Edd = 0.11 for M B h = 1.6 x 10 6 M ; had we chosen the 
radio-loud SED template instead of the radio-quiet one, these 
values would be ~8% higher. Since the SED of UM 625 does 
not, in fact, exactly match the standard shape of the quasar tem- 
plates, an alternative approach is to simply perform a piece - 
wise power-law integration of the observed points (from radio 
to X-ray). This yields a bolometric luminosity that is lower by 
a factor ~5, reducing Lboi/^Edd t0 0.02. Yet a third estimate 
can be obtained from the strength of the Ha emission, which 
scales with the optical continuum emission (e.g., Greene & Ho 
2005b), and hence with Lb i assuming some canonical bolomet- 

11 The specific flux at 5100 A is derived directly from our SDSS spectral fit. If, instead, we estimate ^(5 100 A) from the Ha flux (Greene & Ho 2005b), we obtain 
ckox = -1.76 for the total (narrow plus broad)Hcn flux and a ox = -1.53 for the broad component alone. 
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ric correction for the optical band. Following the formalism of 
Greene & Ho (2007b), the total (broad plus narrow) Ha lumi- 
nosity of 1.0 x 10 41 erg s _1 leads to L bo i = 3.2 x 10 43 erg s _1 
and Lboi/^Edd ~ 0-14. o To summarize, UM 625 is estimated to 
haveL bo i w (0.5-3) x 10 43 erg s" 1 andL bo i/L Edd « 0.02-0.15, 
which is between the values for NGC 4395 and POX 52. 

5.6. Black Hole-Bulge Connection 

We end with a short discussion concerning the connection 
between the mass of the central BH in UM 625 and the proper- 
ties of the bulge of its host galaxy. While the SDSS spectrum 
(Figure 1) does have detected stellar absorption features, we 
are not confident that we have resolved them well enough to 
trust the stellar velocity dispersion derived from the continuum 
fitting (Section 2). In lieu of the stars, we use the velocity dis- 
persion of the narrow-line gas as traced by the low-ionization 
lines [Sn] AA6716, 6731, which are resolved, to estimate the 
stellar velocity dispersion (Greene & Ho 2005a; Ho 2009). For 
a* ~ a[s ii] =68.1 km s _1 , the Mbh-0* relation of Giiltekin et 
al. (2009; their Equation 7) implies M BH = 1.7 x 10 6 M . The 
more recent extension of the Mbh-ct* relation to the low-mass 
end by Xiao et al. (2011) yields M BH = 1.3 x 1O 6 M . Both 
of these estimates are in surprisingly good agreement with our 
virial mass estimate, M B h = 1.6 x 10 6 M , based on the broad 
Ha line (Section 5.1). 

By contrast, and in line with other investigations of low-mass 
AGNs (Greene et al. 2008; Jiang et al. 2011a), the BH in UM 
625 deviates strongly from the BH mass-bulge luminosity rela- 
tion of inactive galaxies. Assuming a spectrum of an Sb galaxy 
(Kinney et al. 1996), the F606W magnitude of the bulge cor- 
responds to My = -19.06 mag, from which the MBH-£y,buige 
relation of Giiltekin et al. (2009; their Equation 8) predicts 
2.2 x 10 7 M , nearly 15 times larger than the virial estimate. 

Now, we know that the stellar population of the pseudobulge 
in UM 625 is quite young, most likely significantly younger 
than the majority of the more massive galaxies used to define 
the MBH-^-buige relation. A proper comparison, therefore, re- 
quires that we apply a correction for age or mass-to-light ratio 
to the bulge luminosity of UM 625. Assuming, as before, an 
Sb galaxy spectrum, for which B—V ~ 0.7 mag (Fukugita et al. 
1995), which, in fact, agrees closely with B — V = 0.67 mag as 
measured by Salzer et al. (1989), a simple stellar population of 
solar metallicity predicts an age of ~ 2 Gyr and a mass-to-light 
ratio M/L v «1.6 M e /L e (Bruzual & Chariot 2003). For a typi- 
cal elliptical galaxy, with B—V ~ 1 .0 mag, M/Ly ^6.3 Mq/Lq . 
Applying this age correction reduces the V-band luminosity by 
about a factor of 4 and the predicted BH mass by about a factor 
of 5. Although the mass discrepancy is now lower, the pseudob- 
ulge of UM 625, like other pseudobulges containing low-mass 
BHs (Greene et al. 2008; Jiang et al. 201 la), systematically de- 
viate from the MeH-^buige relation of classical bulges and ellip- 
tical galaxies. Kormendy et al. (2011) suggest that that there are 
two different modes of accretion: BHs in classical bulges and 
ellipticals grow rapidly via merger-driven gas infall; in contrast, 
BHs hosted in pseudobulges grow mainly via secular evolution 
through slower, stochastic processes. 

Finally, we note that the pseudobulge of UM 625, as in the 
hosts of other low-mass AGNs (e.g. Greene et al. 2008; Jiang 
et al. 2011b), deviates from the Faber- Jackson relation (Faber 
& Jackson 1976) of classical bulges and elliptical galaxies, in 

12 Dong et al. (2012b) give Lboi/^Edd = 0.04 using the same technique, except that they estimate Lboi only using the broad component of Ha, which is 3 times lower 
than the narrow component Ha. Here we use total Ho; while the narrow Ha is seriously contaminated by the star-forming ring as calculated in section 5.3. Our 
strategy is to give a reliable range of Lb i and Lboi/^Edd taking all results into consideration. 



the sense that it has a lower a* for a given luminosity. With 
an /-band absolute magnitude of Mi = -20.27 (assuming, as be- 
fore, a spectrum of an Sb galaxy), the Faber- Jackson relation of 
inactive early-type galaxies (Table 5 of Jiang et al. 2011b) pre- 
dicts a* =124 km s _1 , nearly twice the observed value inferred 
from cr [SII ]. 

6. SUMMARY 

We conducted a comprehensive, multiwavelength study of 
the nuclear and host galaxy properties of UM 625, a type 1 
AGN with a BH mass of 1.6 x 10 6 M determined through the 
detection of broad Ha emission. 

Analysis of Chandra and XMM-Newton observations reveals 
a heavily absorbed (Afa = 9.7^ 8 8 x 10 22 cm~ 2 ) nucleus with an 
intrinsic X-ray luminosity of L2-iokeV = 9.5 x 10 40 erg s _1 . The 
source may be intrinsically X-ray weak (a ox = -1.72) com- 
pared to higher luminosity AGNs. UM 625 belongs to a mi- 
nority of low-mass AGNs detected in the radio, but it is not 
strong enough to qualify as radio-loud. In combination with nu- 
clear photometry at UV, optical, and NIR bands extracted from 
high-resolution HST images, the broad-band SED constrains 
the bolometric luminosity to Lb i ~ (0.5 - 3) x 10 43 erg s _1 and 
W^Edd^ 0.02-0.15. 

We performed a comprehensive analysis of SDSS and HST 
images to quantify the structure and morphology of the host 
galaxy. The galaxy is an isolated, undisturbed, nearly face-on 
SO galaxy with a prominent pseudobulge (Sersic index n = 1 .60) 
that accounts for ~60% of the total light in the R band. The 
pseudobulge has relatively blue colors (B—V ~ 0.7 mag) and is 
mildly disky. Embedded within the central ~ 150-400 pc is a 
UV-bright semi-ring forming stars at a rate of ~0.3 M yr _1 . 
Consistent with other low-mass AGNs, UM 625 follows the 
Mbr-v* relation but not the MeH-^buige relation of inactive 
galaxies. 
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